RET fusions as primary oncogenic drivers and secondary acquired resistance to EGFR tyrosine kinase inhibitors in patients with non-small-cell lung cancer

Background RET fusions are rare oncogenic drivers in non-small cell lung cancer (NSCLC). While activating RET rearrangements are found in NSCLC patients harboring epidermal growth factor receptor (EGFR) genetic alterations at resistance to EGFR inhibitors, the extent to which co-occurring genomic alterations exist and how they might affect prognosis or therapy response is poorly understood. Methods Targeted next-generation sequencing (NGS) was used to assess 380 baseline patients with primary RET fusions and 71 EGFR-mutated NSCLC patients who acquired RET fusions after developing resistance to EGFR-tyrosine kinase inhibitors (EGFR-TKIs). Results Primary RET fusions were more likely associated with females and younger age, with KIF5B being the predominant fusion partner. In baseline patients, both SMAD4 (5.3% vs. 0.0%, P = 0.044) and MYC copy-number gain variants (6.9% vs. 0.0%, P = 0.009) were more frequently co-mutated with KIF5B-RET than CCDC6-RET. By contrast, CDKN2A (11.3% vs. 2.4%, P = 0.003) mutations were significantly enriched in CCDC6-RET-rearranged baseline patients. A significant increase in the proportion of CCDC6-RET was observed in acquired RET-rearranged patients (47.3% vs. 22.5%, P < 0.001). The median progression-free survival (PFS) of patients harboring RB1 and TP53 double-mutations (5.5 vs. 10.0 months, P = 0.020) or ERBB2 amplification (5.6 vs. 10.0 months, P = 0.041) was significantly shorter than the wild-type counterparts. Moreover, we identified that RET fusions were more likely associated with acquired resistance (AR) to third-generation EGFR-TKIs than previous generations of EGFR-TKIs. Conclusions In conclusion, we depicted the mutational profiles of NSCLC patients who harbor RET fusions at baseline or after resistance to EGFR-TKIs. Furthermore, our results suggest that RET fusions mediate secondary resistance to third-generation EGFR-TKIs and might be associated with poor prognosis in patients with NSCLC. Supplementary Information The online version contains supplementary material available at 10.1186/s12967-022-03593-3.

The fusion of RET with another unrelated gene occurs due to an aberrant DNA repair process [2]. The resulting fusion product activates various downstream signaling pathways that play essential roles in cell proliferation and survival [5]. Previous studies have shown that RET fusion-positive (RET+) NSCLC patients testing negative for EGFR/ALK/BRAF/ROS1 are usually young never-smokers with ADC [6][7][8][9]. While the most common RET fusion partners in NSCLC patients are KIF5B and CCDC6, other reported partners include NCOA4, TRIM33, etc. [10]. The overall survival of CCDC6-RET+ baseline patients was nearly three times longer than those with KIF5B-RET fusions (median: 113.5 vs. 37.7 months, P = 0.009) [11]. Besides, treatment responses of RET inhibitors were heterogeneous among baseline RET+ patients harboring different fusion variants [9,[11][12][13][14], highlighting the importance of fusion partner types in clinical outcomes toward targeted therapy.
More recently, receptor tyrosine kinase (RTK) fusions have emerged as a rare but targetable acquired resistance (AR) mechanism in EGFR-mutated NSCLC patients on EGFR-TKI treatment. Notably, RET fusions are the most commonly reported RTK fusions that mediate AR to EGFR-TKIs [15]. Within secondary RET fusions, CCDC6-RET is the most common fusion variant, followed by NCOA4-RET. Interestingly, NSCLC patients harboring KIF5B-RET fusions showed minimal response after RET TKI (RXDX-105) treatment, whereas the response rate was 67% in non-KIF5B-RET+ NSCLC patients [16]. Dual blockade of EGFR driver mutation and RET fusion, such as CCDC6-RET and NCOA4-RET, demonstrated safety and clinical efficacy in both clinical and preclinical studies [17].
To date, two highly potent RET-specific TKIs, selpercatinib and pralsetinib, have been approved by the US Food and Drug Administration (FDA) for the treatment of advanced or metastatic RET-altered NSCLC and thyroid cancers. Selpercatinib and pralsetinib effectively against RET alterations, including CCDC6-RET and KIF5B-RET fusions, RET activating mutations (C634W and M918T), and RET gatekeeper mutations V804L/M/E [18,19]. Remarkably, selpercatinib has > 100-times selectivity against VEGFR2, and pralsetinib has 87-times selectivity against VEGFR2 and 20-times selectivity against JAK1 [20]. Findings from the phase I/ II LIBRETTO-001 trial (NCT03157128) demonstrated that selpercatinib has an overall response rate (ORR) of 64% in previously treated NSCLC patients and ORR of 85% in treatment-naïve RET-altered NSCLC patients [21]. In addition, the antitumor potential of selpercatinib is irrespective of specific RET fusion types. On the other hand, initial data from the phase I/II ARROW trial (NCT03037385) demonstrated that pralsetinib has a high potency and durable activity and is well-tolerated in adult patients with metastatic RET-altered NSCLC. The ORR in previously treated patients was 61%, and ORR in treatment-naïve patients was 70% [22]. Most treatment-related adverse events (TRAE) of selpercatinib and pralsetinib are mild and controllable, including anemia, elevated alanine aminotransferase and hypertension [5,20]. However, the safety profiles of these two RET inhibitors need to be further studied, given that some safety warnings have been reported [20][21][22][23][24].
In this study, we delineated the mutational profiles of 380 baseline and 71 EGFR-mutated NSCLC patients who acquired RET fusions after resistance to EGFR-TKIs by targeted NGS and revealed RET fusion partners associated with primary and acquired patients. We also investigated the impact of co-occurring genetic alterations in RET-rearranged NSCLC patients, which might explain the poor prognosis of patients harboring secondary RET fusions.

Patients and sample collection
Tumor tissue and/or plasma samples were collected from 451 RET+ NSCLC patients admitted to all participating hospitals between June 2015 and June 2021. Specifically, formalin-fixed paraffin-embedded (FFPE) tumor or fresh tumor tissue were confirmed by pathologists from the centralized clinical testing center. 5-10 mL of peripheral blood was collected from each patient in EDTA-coated tubes (BD Biosciences) and shipped to the clinical testing center within 48 h of blood collection for the following tests. Clinical characteristics and treatment history were extracted from medical records. Progression-free survival (PFS) was defined as the time from the initiation of the treatment to disease progression/patient death. Patients who had not progressed were censored at the date of their last follow-up. This study was conducted in accordance with the declaration of Helsinki and was approved by the Ethical Review Board of Fujian Medical University Union Hospital. Informed written consent was obtained from each subject before sample collection.

Targeted next-generation sequencing
DNA extraction, library construction, and targeted NGS were performed as previously described in a Clinical Laboratory Improvement Amendments (CLIA)-certified and College of American Pathologists (CAP)-accredited clinical testing laboratory (Nanjing Geneseeq Technology Inc., Nanjing, China) [25,26]. FFPE samples were de-paraffinized with xylene followed by genomic DNA extraction using QIAamp DNA FFPE Tissue Kit (Qiagen Cat. No. 56404) according to the manufacturer's instructions. Genomic DNA from fresh tumor tissue was extracted using the DNeasy Blood and Tissue Kit (Qiagen Cat. No. 69504) according to standard protocols. Peripheral blood samples were centrifuged at 1800g for 10 min, followed by cell-free DNA (cfDNA) extraction and purification using QIAamp Circulating Nucleic Acid Kit (Qiagen Cat. No. 55114). Genomic DNA of white blood cells in sediments was extracted using the DNeasy Blood and Tissue Kit (Qiagen Cat. No. 69504) as normal control. Genomic DNA was qualified using Nanodrop2000 (Thermo Fisher Scientific, Waltham, MA), and cfDNA fragment distribution was analyzed on a Bioanalyzer 2100 using the High Sensitivity DNA Kit (Agilent Technologies, Santa Clara, CA, 5067-4626). DNA quantification was performed using the dsDNA HS assay kit on a Qubit 3.0 fluorometer (Life Technology, US). NGS libraries were prepared using the KAPA Hyper Prep kit (KAPA Biosystems) with an optimized manufacturer's protocol for different sample types. Targeted capture enrichment was performed as previously described [27]. The target-enriched library was then sequenced on HiSeq4000 or HiSeq4000 NGS platforms (Illumina) according to the manufacturer's instructions.

Mutation calling
Sequencing data were processed as previously described [25]. In brief, the data was first demultiplexed and subjected to FASTQ file quality control using Trimmomatic [28]. Low-quality data (QC below 15) and N bases were removed. Raw reads were then mapped to the Human Genome (hg19) using Burrows-Wheeler Aligner (BWA-mem, v0.7.12; https:// github. com/ lh3/ bwa/ tree/ master/ bwakit). Genome Analysis Toolkit (GETK 3.4.0; https:// softw are. broad insti tute. org/ gatk/) was employed to perform local realignment, base quality score recalibration, and detect germline mutations. Picard was used to remove PCR duplicates. VarScan2 was applied to detect single-nucleotide variations (SNVs) and insertion/deletion mutations. The limit of detection (LOD) of tumor tissues and plasma samples under specific sequencing depths has been tested repetitively in Nanjing Geneseeq Technology Inc. to ensure the mutation calling results are consistent and have optimal performance for identifying genetic alterations. According to its internal specifications, sequencing of tissue and plasma ctDNA via targeted NGS can reach a sensitivity of 98% and a positive predictive value (PPV) of 95% [26,29,30]. SNVs were filtered out if the VAF was less than 1% for tumor tissue and 0.3% for plasma samples. Common SNVs were excluded if they were present in > 1% population in the 1000 Genomes Project or the Exome Aggregation Consortium (ExAC) 65,000 exomes database. The resulting mutation list was further filtered by an in-house list of recurrent artifacts based on a normal pool of whole blood samples. Parallel sequencing of matched white blood cells from each patient was performed to remove sequencing artifacts, germline variants, and clonal hematopoiesis. Genomic fusions were identified by FACTERA [31] with default parameters (≥ 2 reads). The fusion reads were manually reviewed and confirmed on Integrative Genomics Viewer (IGV). By definition, RET fusions were annotated as a protein fusion involving: (i) a 5′ non-RET partner gene and (ii) an intact 3′ RET kinase domain (NM_020975: exon [12][13][14][15][16][17][18]. In this study, RET fusions were classified into canonical (single or compound KIF5B-RET and CCDC6-RET) and noncanonical RET. In this regard, a noncanonical RET fusion includes: (i) a rearrangement with a partner rather than KIF5B and CCDC6; (ii) A rearrangement with a novel partner gene; and (iii) a rearrangement with an intergenic region (IGR). Tumor mutational burden (TMB, mutation per Megabase) was determined based on the number of missense mutations in the targeted regions of the gene panel covering 0.85 Mb of coding genome, excluding known driver mutations as they are overrepresented in the panel. Chromosome instability score (CIS) was defined as the proportion of the genome with aberrant (purity-adjusted segment-level copy number ≥ 3 or ≤ 1) segmented copy number [32]. For gene-level analyses in the baseline cohort, genetic alterations with a mutation frequency ≥ 2% were considered frequently mutated. All 139 genes in the PULMOCAN ™ gene panel (Geneseeq Technology Inc.) were subjected to pathway enrichment analyses. Specifically, the genetic alterations identified by targeted panel sequencing were first categorized into ten signaling pathways associated with common hallmarks of cancer that control cell-cycle progression, apoptosis, cell proliferation and growth [33]. The proportion of baseline patients with specific RET fusion variants harboring mutation(s) in the relevant pathways was compared to reveal the differences of co-existing mutations in baseline RET+ patients.

Statistical analysis
Fisher's exact tests were used to compare the categorical variables between groups. Kaplan-Meier curves were used to analyze the PFS of various patient groups, and the statistical difference was analyzed using the log-rank test. A two-sided P value of less than 0.05 was considered significant for all tests unless indicated otherwise (*P < 0.05, 0.01 < **P < 0.05, ***P < 0.001). All statistical analyses were performed using R (version 4.1.2).

Distribution of RET fusion partners between primary and acquired RET+ patients
By targeted NGS, a total of 491 RET fusions were identified in 451 RET+ patients (Fig. 2a), of which KIF5B-RET (51.1%) was the most frequently observed RET rearrangement, followed by CCDC6-RET (23.4%). Notably, 39 patients (8.6%) carried more than one RET fusion, of whom 36 were baseline and 3 were acquired RET+ patients. Next, we tried to depict functional domains and breakpoints of canonical RET fusions. In particular, 78%  were EGFR-mutated patients who acquired RET fusions after resistance to EGFR-TKIs. Only baseline RET fusion-positive (RET+) patients with available FFPE and/or biopsy tumor samples were qualified for the following tests. Acquired RET+ patients (N = 71) were those who gained RET fusions after EGFR-TKI treatments targeting primary EGFR oncogenic mutations. Samples belonging to this category were divided into three groups depending on when RET fusion was detected. Specifically, group 1 contains 13 patients treated with first-line (1L) 1st-or 2nd-G EGFR-TKIs. Group 2 includes 51 patients who received second-line (2L) 3rd-G EGFR-TKI treatment. Lastly, group 3 consists of 7 patients previously treated with 1L 3rd-G EGFR-TKIs within intron 11 of RET, while translocation events may also occasionally occur within intron 7 and 10, resulting in the inclusion of the RET transmembrane domain [34]. The next question we wanted to address was whether the percentage of specific RET fusion variants differs in baseline and acquired RET+ patients. We identified 417 and 74 RET rearrangements in these two cohorts, respectively (Fig. 2c, d). The most commonly observed fusion variants in baseline patients were KIF5B-RET (53.7%), followed by CCDC6-RET (22.5%, Fig. 2c). Surprisingly, KIF5B-RET accounted for only 8.1% of total translocation events in baseline EGFR-mutated NSCLC patients resistant to EGFR-TKIs, whereas CCDC6-RET was identified in 47.3% of translocations events, ranking as the predominant RET fusion variant in acquired RET+ patients (Fig. 2d). The identification of CCDC6-RET was the most common fusion variant, followed by NCOA4-RET, was consistent with previous findings [15]. Furthermore, these results implied that KIF5B-RET and non-KIF5B-RET fusions might have different functionalities in NSCLC. At the same time, noncanonical RET fusions, including those with hitherto unreported partner genes [10], were identified in 23.8% and 44.6% of baseline and secondary patient cohorts, respectively (Additional file 1: Table S1).
As RET fusion variants distributed differentially in baseline and acquired RET+ patients, we investigated whether RET fusion types could be associated with the patient's clinical manifestations. Interestingly, females were more likely associated with KIF5B-RET (P = 0.006) and noncanonical RET fusions (P = 0.015) than CCDC6-RET in baseline patients (Additional file 2: Fig. S1a). However, neither the patient's age nor cancer stage at diagnosis was directly associated with RET fusion types in baseline or secondary patients (Additional file 2: Fig.  S1b-e).  (Fig. 3a).

RET fusion confers a secondary resistance mechanism to EGFR-TKIs
Despite their rarity, it is clear from previous studies that RTK fusions, such as RET rearrangements, are actionable resistance mechanisms to EGFR-TKIs. We aim to increase awareness of this emerging paradigm by comprehensively profiling baseline EGFR-mutated NSCLC patients who acquired RET fusions after developing resistance to EGFR-TKIs. In acquired RET+ patients, second-site EGFR mutations, such as T790M (50.7%), C797S/G (16.9%) and L718V/Q (5.6%), were among the top frequently mutated gene alterations (Fig. 4a). It was also apparent that cell cycle control pathway gene alterations, such as TP53, RB1, CDKN2A, CDKN2B, CDKN1B, and CDK6, co-existed with acquired RET fusions. In addition, ALK and NTRK1 fusions that belong to the RTK family were identified to co-occur with RET fusion in two patients. These results implied that RTK fusions and other genetic variants could potentially serve as resistance mechanisms to anti-EGFR treatment in NSCLC.
To characterize the AR mechanism through RET fusions in EGFR-mutated NSCLC patients and their survival outcomes, we compared the PFS among patients treated with different EGFR-TKI regimens. However, no significant difference was observed (Additional file 5: Fig. S4a). Alternatively, as 2L 3rd-G EGFR-TKI treatment has been employed as a standard practice for EGFR-mutated NSCLC patients with progressive disease on first-line targeted therapy, we performed integrated survival analyses using group 2 patients. Notably, we excluded one patient with the rare EGFR double-mutants and one with SCC histology, leaving 49 patients with only EGFR 19-Del or L858R mutations in the refined cohort (Additional file 5: Fig. S4b). We first performed a univariate analysis on various demographic and mutational features of refined cohort 2 patients. Interestingly, bypass pathway genetic alterations (KRAS and PIK3CA activation mutations, ERBB2 and MET amplification, ALK and NTRK fusions), RB1 and TP53 co-mutation, and ERBB2 copy-number gain had significant associations with the prognosis of acquired RET+ patients who underwent 2L 3 rd -G EGFR-TKI treatment (Additional file 1: Table S2). We then performed multivariate analysis using these three features, and we discovered that they were no longer statistically significant, with bypass mutations and RB1and TP53 co-mutations being close to significance. This result suggested that there might be some interactions or correlations among these prognostic features. Nonetheless, the PFS of patients harboring bypass pathway gene alterations was significantly shorter than wild-type patients (median: 6.7 vs. 12.0 months, P = 0.020, Fig. 4b). In addition, the PFS of patients with TP53 and RB1 co-mutations was significantly shorter than wild-type counterparts (median: 5.5 vs. 10.0 months, P = 0.02, Fig. 4c). ERBB2 gene amplifications that were present in 5.6% of acquired RET+ patients were also associated with poor prognosis in NSCLC patients (median: 5.6 vs. 10.0 months, P = 0.041, Fig. 4d). Notably, although the majority of acquired RET fusions were detected by circulating nucleic acid method using plasma ctDNA, 4 out of 49 patients in the refined cohort contributed only FFPE samples, raising a concern of missing rare fusions compared to tumor-based comprehensive genomic profiling. We, thereby, excluded patients with only tissue samples and re-performed the analysis in Fig. 4 using only ctDNA data (Additional file 6: Fig. S5a). Overall, our results about the co-occurring genetic alterations were consistent when using either all samples (tissue and ctDNA) or only ctDNA (Additional file 6: Fig. S5bd), suggesting that our conclusions are not likely to be affected by different sample types. Collectively, comprehensive molecular profiling of acquired RET + patients showed that co-existing genomic alterations, such as TP53 and RB1 co-mutations and ERBB2 amplification, might be associated with a poor prognosis in EGFR-mutated NSCLC patients upon drug resistance.

RET fusions are more frequently associated with third-generation EGFR-TKIs
Previous case reports have linked RET fusions, such as CCDC6-RET, TRIM24-RET, NCOA4-RET and ERC1-RET, to osimertinib resistance in EGFR-mutated NSCLC patients [35][36][37][38][39]. Here, we demonstrated a higher incidence rate of RET fusions in patients who developed AR to 3rd-G EGFR-TKIs from a population perspective. Interestingly, the incidence of RET fusions was highest in patients who underwent 2L 3rd-G EGFR-TKI treatment (1.5%, 51/3330), followed by patients treated with 1L 3rd-G EGFR-TKIs (1.2%, 7/589, Fig. 4e). These numbers were about ten-fold higher than those in patients treated with first-or second-generation EGFR-TKIs at the front line (0.15%, 13/8732). Overall, we demonstrated that RET fusions were more likely associated with EGFR-mutant NSCLC patients who received therapeutic interventions targeting EGFR with third-generation EGFR-TKIs. The complexity of non-RET secondary mutations in acquired RET+ patients might contribute to the differential treatment responses to follow-up therapies.

Discussion
With comprehensive genomic profiling, clinicians can use the knowledge of specific clinical features associated with individual alterations to optimize therapeutic decision-making. Here, we retrospectively analyzed the mutational profiles of 451 patients carrying either baseline or acquired RET fusions to characterize the roles of RET fusions in NSCLC.
As a rare oncogenic driver mutation, RET rearrangement occurs in 1-2% of NSCLC patients [1,2]. It has been shown that RET-rearranged NSCLC patients are more likely associated with ADC. Consistent with this notion, we found that 76.1% of patients developed ADC as the predominant histological type in the baseline cohort. On the other hand, a much-debated question is whether baseline RET fusions can correlate with the patient's gender or age [8,40,41]. Hence, we compared clinical characteristics of baseline patients with RET fusions, such as gender and age. Our results showed that primary RET fusions were more likely to occur in females than males (55.3% vs. 44.7%, P = 0.038) and were associated with younger patients (58.9% vs. 41.1%, P < 0.001). The controversy could be due to variations in the subject's ethnicity, genetic background, environmental factors, and even lifestyles. Then, we compared whether the distribution of specific RET fusions differed in baseline and acquired RET+ patients. We found that the predominant fusion type in baseline patients was KIF5B-RET, whereas CCDC6-RET was most frequently identified in patients who acquired RET fusions at resistance to EGFR-TKIs. The proportion of patients harboring NCOA4-RET fusions also increased from 1.0 to 16.2% in acquired RET+ patients. These results were consistent with previous findings, suggesting that KIF5B-RET and non-KIF5B-RET fusions might have different functionalities in EGFR-TKI progression [17,35]. Lastly, we characterized mutational profiles of the 380 baseline patients and demonstrated that concurrent gene alterations, such as SMAD4 mutations and MYC copy-number gain, were more frequently associated with KIF5B-RET than CCDC6-RET fusions in baseline RET-rearranged NSCLC patients. Our results may provide insights into why NSCLC patients harboring KIF5B-RET fusions have a nearly three times shorter overall survival than those harboring CCDC6-RET fusions.
The diversity and complexity of molecular mechanisms underlying the acquired adaptation of cancer cells to targeted therapies, such as EGFR-TKIs, is an area of active investigation. In this study, we demonstrated that RET fusions, as a rare but actionable AR mechanism to EGFR-TKIs, confer a poor prognosis in EGFR-mutated NSCLC patients. It has been previously suggested that patients harboring TP53 and RB1 co-mutations are at unique risk of histologic transition from ADC to SCC or eventually small cell transformation [42]. Here, our results showed that RET+ patients with co-occurring TP53 and RB1 double-mutations had significantly shorter PFS than wild-type patients, highlighting the role of TP53 and RB1 in controlling cell proliferation and disease progression. In addition, ERBB2 copy-number gain can also present a similar effect, resulting in reduced PFS in patients who acquired RET fusions at resistance to EGFR-TKIs. It should also be noted that prognostic-related factors being examined in the survival analyses, including bypass pathway gene alterations, TP53 and RB1 co-mutations, and ERBB2 copy-number gain, were not completely independent. The limited sample size may have impaired the statistical power to achieve significant results. Further studies using larger cohorts that consider additional variables need to be undertaken to better understand prognostic factors associated with the patient's survival.
In the final part of our study, we compared the incidence of secondary RET fusions in NSCLC patients who progressed on different EGFR-TKI therapies. A ten-fold higher incidence of RET fusions was observed in patients who underwent third-generation EGFR-TKI treatment than those treated with front-line 1st-/2nd-G EGFR-TKIs. This result was consistent with a previous finding which suggested that RET fusions were significantly enriched after 3rd-G EGFR-TKI treatment [15]. As previously mentioned, selpercatinib and pralsetinib have been granted FDA approval, showing equipotent for treating RET-rearranged NSCLC and thyroid cancer with minor and controllable adverse effects. The ORR of previously treated RET+ NSCLC patients on selpercatinib and pralsetinib were 64% and 61%, respectively, while the median PFS of previously treated patients on these two RET inhibitors was 16.5 months and 17.1 months, respectively [21,22]. Given the high efficacy and mild side effects of the two RET-specific inhibitors, it is worth investigating their clinical utility to treat EGFR-mutated NSCLC patients who acquired RET oncogenic alterations after TKI resistance, especially those receiving 3rd-G EGFR-TKIs. Despite efforts we made to systemically characterize the mutational profiles of baseline and acquired RET+ patients, there are limitations in our study. We did not validate the authenticity of RET fusions with hitherto unreported partner genes. However, we reasoned that patients without known driver mutations might harbor functional RET fusions since RET fusions usually occur mutually exclusively. Due to the unavailability of RET-specific inhibitors during the patient's treatment, we could not evaluate the PFS of baseline patients with different RET fusion partners on RET-specific inhibitors or patients with acquired RET fusions after EGFR-TKI resistance. However, dynamic monitoring of these acquired RET+ patients is currently ongoing. We intend to investigate the treatment outcomes for EGFR-mutated NSCLC patients harboring acquired RET fusions on the follow-up RET-specific inhibitor therapy upon obtaining more data on these patients.

Conclusions
In conclusion, we systematically evaluated the mutational profiles of RET+ baseline patients and those who acquired RET fusions at secondary resistance to EGFR-TKIs. We identified unique genetic features explaining the differential treatment responses in NSCLC patients harboring baseline RET fusions. Furthermore, we demonstrated that 3rd-G EGFR-TKIs are more likely associated with secondary RET fusions. The high efficacy and mild side effects of the two RET-specific inhibitors may provide treatment options for EGFR-mutated NSCLC patients who developed RET fusions following TKI resistance, especially those on 3rd-G EGFR-TKI treatments.